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Abstract 

We consider the corrections to deuterium hyperfine structure originating from the 
two-photon exchange between electron and deuteron, with the deuteron excitations in the 
intermediate states. In particular, the motion of the two intermediate nucleons as a whole 
is taken into account. The problem is solved in the zero-range approximation. The result 
is in good agreement with the experimental value of the deuterium hyperfine splitting. 
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1 Introduction 



The hyperfine (hf) splitting in deuterium ground state has been measured with high accuracy. 
The most precise experimental result for it was obtained with an atomic deuterium maser and 
constitutes pQ 

i/ ex p = 327 384.352 522 2(17) kHz. (1) 
Meanwhile the theoretical calculation including higher order pure QED corrections gives 

z/ Q ed = 327 339.27(7) kHz. (2) 

The last number was obtained by using the theoretical result for the hydrogen hf splitting from 
Ref. 

1420 451.95(14) kHz 

which does not include proton structure and recoil radiative correction, and combining it with 
the theoretical ratio of the hf constants in hydrogen and deuterium from Ref. [B] 

4.339 3876(8) 

based on the ratio of the nuclear magnetic moments and including the reduced mass effect 
in \ip(0)\ 2 . 

It was recognized long ago that the discrepancy 

^exp - ^qed = 45 kHz (3) 

is due to the effects caused by the finite size of deuteron. Such effects are obviously much larger 
in deuterium than in hydrogen. Corresponding contributions to the deuterium hf splitting were 
discussed long ago with some intuitive arguments |1J, and then in more detail in Refs. |S1IS1IZ|- 

We believe that in the past the most systematic treatment of such effects, which are due to 
the electron-deuteron interaction of second order in a, was performed in Ref. The effective 
hamiltonian of the hf interaction of second order in a = e 2 /4ir was derived therein from the 
elastic forward scattering amplitude of virtual photons off the deuteron. 

In particular, the low-energy theorem for forward Compton scattering [01 E3 EJ Ej was 
generalized in jH] to the case of virtual photons and a target with arbitrary spin. The cor- 
responding contribution of the momentum transfers k, bounded from above by the inverse 
deuteron size x = 45.7 MeV, to the relative correction to the deuterium hyperfine structure is 

a ,i 3a / 3 \ m e , x 

K= 7r(^- 2 —In—. (4) 

87r V /W rn p m e 

Here m e and m p are the electron and proton masses, respectively, = 0.857 is the deuteron 
magnetic moment. The relative corrections A's are defined here and below as the ratios of the 
corresponding contributions to the ed scattering amplitude to the spin-dependent Born term 
in this amplitude, 

Itxol 

To = M°" ■ s )> (5) 
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where s is the deuteron spin. 

At larger momentum transfers, k > x, the amplitude of the Compton scattering on a 
deuteron is just the coherent sum of those amplitudes on free proton and neutron. This cor- 
rection equals 

A£= ^-(^-2^-3 + ^)^ln^. (6) 

4-7T Hd m p x 



Here fi p = 2.79 and \L n = — 1.91 are the proton and neutron magnetic moments; 
m p = 770 MeV is the usual hadronic scale. 

Strong numerical cancellation between A e i and A in is worth mentioning. 

The next correction to the deuterium hyperfine structure (hfs), obtained in Ref. [8], is 
induced by the deuteron virtual excitations due to spin currents only. It is 

. (!) _ 3a Op - /i TO ) 2 me, mp (r? s 

in o ' \ / 
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There is also a correction due to a finite distribution of the deuteron charge and magnetic 
moment 3 . In the zero-range approximation used in Ref. [S] this correction is 

A f = -a-^(l + 21n2). (8) 

O X 

2 Leading inelastic nuclear correction 
to deuterium hyperfine structure 

Leading inelastic nuclear correction is on the relative order a m e / x (as well as the Zemach 
correction (JBJ). The corresponding effect calculated in Ref. |H] is additionally enhanced by a 
large factor \i v — fi n = 4.7. In the present paper we consider two more effects on the same order 
am e j k. Though both of them are proportional to \x v + fi n = 0.88 (and thus are essentially 
smaller numerically than that considered in [S|), we believe that their investigation is worth 
attention. 

We use the gauge A = where the photon propagator is 

dim kik m 

Di m (uJ, k) = — zr , di m — Sim — — 5~ ] Dqq = Do m = 0. (9) 

uj 1 — k 2 w 

The electron-deuteron nuclear-spin-dependent scattering amplitude generated by the two- 
photon exchange is 

/d 4 k dj m dj„ 7,-(Z — k + rn P )7,- 
(2^ fc 4 k 2 - 2ik (10) 

Here 1^ = (m e , 0, 0, 0) is the electron momentum. The structure t,-(Z — k + m e )^j reduces to 
— iuj eijiai where a is the electron spin. We calculate the nuclear matrix elements entering the 
deuteron Compton amplitude M mn in the zero-range approximation (zra) which allows us to 
obtain all the results in a closed analytical form. 

The inelastic 1/x contribution to hfs is induced by the combined action of the convection 
and spin currents. Since the convection current is spin-independent, all the intermediate states 
are triplet ones, as well as the ground state. Therefore, here the spin current operator 

g 

ik x [fipcr p exp(ikr p ) + fi n er n exp(zkr n )] 



2m p 



simplifies to 

g 

z[k x s] [/i p exp(?kr p ) + //„exp(zkr„)]. (11) 



2m p 



3 In the case of hydrogen this problem was considered many years ago by Zemach |13j . 



In the initial state | 0) the deuteron is at rest. But in the excited state the system of nucleons as 
a whole moves with the momentum k, so that its wave function is | n) exp(zkR), where | n) refers 
to the deuteron internal degrees of freedom and is a function of r = r p — r n ; R = (r p + r n )/2 
is the deuteron centre of mass coordinate. Thus, a typical matrix element of the spin current 
can be written as 

g 

ik x (n| exp (— ikR)s [/i p exp(ikr p ) + jj n exp(ikr n 



2m p 



— - — ik x (n|s[yU p exp(zkr/2) + /i n exp(— zkr/2)]|0) . (12) 

As to a typical matrix element of the convection current, it transforms now as follows: 

g 

- — (n| exp(— ikR)p B exp(ikr p ) + exp(ikr„)pJ0) 

2m p ' ' 

(n|(p+ -)exp(ikr/2) + exp(zkr/2)p|0) = — (n|pexp(ikr/2)|0) ; (13) 



2m p 2 ' ' m p 

here p acts on the relative coordinate r. 

At first let us take as intermediate states \n) in the corresponding nuclear Compton am- 
plitude just plane waves, eigenstates of p. In this way we take into account all the states 
with / ^ 0, which are free ones in our zero range approximation, and in addition the 3 Si wave 
function in the free form ijj p (r) = smpr/pr (the deviation of the 3 Si wave function from the 
free one will be considered below). Then, with the zra deuteron wave function 

. , . f~x exp(-xr) . . 

V 27r r 

the only matrix element entering the amplitude is 



^|exp( ±i kr/2)|p) = (p± gr +x2 . (15) 



In this way the amplitude itself simplifies to 



2 
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dp ( fi p jj n 
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2m p J J ir 2 U(p-k/2) 2 + x 2 ] 2 [(p-k/2) 2 + x 2 ][(p + k/2) 2 + x ; 

2p m i € nrs k r s s 2p n i € mrs k r s s mc\ 
X uj 2 - (p 2 + P/4 + x 2 ) 2 /m 2 ' 1 j 

Due to the account for the motion of the system as a whole in the intermediate states, 
this expression differs from the corresponding one from our previous paper jH] in two respects. 
First, in jHj the operator p p in (|T3|) was identified with p. Thus, therein instead of 2p m ^ n in the 
analogue of the present formula (|TK|). we obtained {2p — fc/2) m>n . At present, in (fTH|) the term 
proportional to /i n is an odd function of p, and therefore vanishes after integration over dp. 
Second, in the denominator the energy difference has acquired the contribution k 2 /4m p , which 
is the kinetic energy of the pn system as a whole, and thus (p 2 + x 2 )/m p has transformed into 
(p 2 + k 2 /4 + x 2 )/m p . 

Now we substitute (|16|) into (JIUj) and take the integral over uj under the condition uj ^> h 2 jm. 
For the relative correction to the hf structure we obtain 



( 2 ) 2ax/i p m e f f dp dk pk 



■m 



7r^ d m p JJ [(p - k/2) 2 + x 2 ] 2 



nip 3 



p 2 + fc 2 /4 + x 2 2k 



(17) 



The result of integration over p and then over k reads 



2 p p m e 6a fi p m e m p 

A|„ = a — — In — . (18) 

Hd X 7T fj, d m p X 

The logarithmic contribution here originates from integration of the term 3 /2k in the square 
brackets in (|17|) over the range x 2 /m p « C x. The result (|18|) differs from the corresponding 
one of jS] by a term proportional to p p + p n , which is relatively small numerically. It is only 
natural that our present account for the motion of the proton-neutron system as a whole in the 
intermediate states results in the correction proportional to fi p + p n . 

Let us calculate now the correction A\ n ' corresponding to the effect of deviation of the 
intermediate 3 Si wave function *& p (r) from the free one. In the zra ^ p {r) reads 



sinpr 1 exp(ipr) x sin pr — p cos pr 

^pv) = ; : = 7 ; r~v • (l^J 

pr x + ip r pr{x + %p) 



It follows, for instance, from the orthogonality to the deuteron wave function ()14j) . We use 
below the function 

i \ t / \ t * / \ , / \ / * / \ pcosp(ri + r 2 ) - xsinp(ri + r 2 ) , nn . 

Pp (ri,r 2 ) = y p ( ri )%(r 2 ) - ^(n)^(r 2 ) = (^+^W 2 ' (20) 

(3) 

Then, after integration over uj the expression for A in reads 
( 3 ) _ 4a + /i„)m, 



A L = 3 / / dkdpp I I dT 1 dT 2 ipo(ri)ipo(r2)pp(r 1 ,r 2 ) 

7T HtjTflp 



sin kr 2 (1 + xr 2 ) f sin kr 2 
kr 2 (kr 2 ) 2 \ kr 2 



m p 3 



p 2 + k 2 /A + x 2 2k 



(21) 



sin kv\ 
kr\ 

The integral over p is 

dp p^A + l 2 = ^k~ 2 [(g + X) eXp[ " g(ri + r2)] - + ' ( 22 ) 

where Q = a/x 2 + /c 2 /4. Now we integrate (|21|) over r 1; r 2 , and then over fc. The final result 
for the discussed correction reads: 

A (3) = -q /ip + /i "^l(2-21n2)+ ^ /ip + /i "^ln^. (23) 
/id x 3 7r p d m p x 

Again, it is only natural that due to common selection rules the contribution of the 3 Si interme- 
diate state is proportional to /i p + /i n . Let us note also that the first term in (|2~3j) is additionally 
suppressed by a small numerical factor (2 — 21n2)/3 = 0.20. 

At last, we wish to come back to the effect due to a finite distribution of the deuteron charge 
and magnetic moment. The Zemach correction Af (jHJ can be also easily derived in the present 
approach. Using the identity 

(V>o |pexp(ikr/2)|^o) = ^(^o I exp(ikr/2)|^o) , 
we obtain the corresponding amplitude: 

ML = bp + K»)"l F W- 1 ! ^- (p2 + fcV 4 + x 2 ) 2 /m 2 ' (24) 



where F(k) is the deuteron form factor in the zero-range approximation 



F(k) — (ipo I exp(ik.r/2)\il> ) — — arctan — . (25) 

k 4x 

Let us note that in our approximation both electric and magnetic form factors, which in the 
present case enter the convection current and spin current matrix elements respectively, coincide 
and equal F{k). 

The integration over u leads to the following result for the relative correction Af 

8a(p n + u„)m P f°° dk r _ , , u„ + u„ m e 1 . , . 

A f — — \tlZJ^l / — \F (k) — 1 = -a P — - l + 21n2 . 26 

There is no logarithmic term in Ay since F 2 (k) — 1 ~ k 2 for « x. In fact, the result (|2l)j) 
agrees with (jSJ) since to our accuracy fi p + fi n = fi d - 

The corrections (|18|). ()23|). combine into a compact result 

A c = - a — — ^ — — — In . 27 

fjL d k ix fx d m p x 

Let us notice that its logarithmic part coincides with the corresponding logarithmic term in 
|H] (see formula (27) therein). This is quite natural: the log is dominated by small k, so that 
extra power of k, which arises from the recoil of the pn system as a whole, cannot influence it. 

The leading term in (j2*7|) coincides with the result of Ref. [3] 4 . However, we could not find 
any correspondence between our arguments and those of Ref. [Hj- In particular, it is stated 
explicitly in Ref. [5 s that the motion of the intermediate pn system as a whole is neglected 
therein. 



3 Discussion of results 

Our total result for the nuclear-structure corrections to the deuterium hf structure, comprising 
all the contributions, (gj), ©, ©, (j2TJ) . is 

fi n m e 3a fi p - fi n m e m p 
A = — a — In — - (28) 

£*<f X 7T fl d m p X 

3a (fip - fi n ) 2 m e ln ^»p 
87r fid m p x 

3a 1 . , . m e , x 3a 1 . 9 , . m e , m 

+ T — (t4 - - 3 — In— + — — (fj 2 - 2fx p -3 + nD-t^-t. 

87T pd Kip ™>e 47T /irf v VTLp X 

Numerically this correction to the hf splitting in deuterium constitutes 

Av = 50 kHz. (29) 

It should be compared with the lacking 45 kHz (see 0)- We believe that the agreement 
is quite satisfactory if one recalls the crude nuclear model (zra) used here; in particular the 
deuteron form factors, as calculated in the zra, are certainly harder than the real ones, and 
thus in zra the negative Zemach correction is underestimated. 

4 We are sorry for misquoting the result of Ref. |S] in our paper [S]. 



Let us mention here that in a recent paper ^3] elastic contributions and the Zemach effect 
were considered in a quite different theoretical technique, but with some phenomenological 
desription of the deuteron form factors. The result is smaller than the corresponding part of 
ours by 13 kHz. 

Clearly, the nuclear effects discussed are responsible for the bulk of the difference between 
the pure QED calculations and the experimental value of the deuterium hf splitting. The 
calculation of this correction, including accurate treatment of nuclear effects, would serve as 
one more sensitive check of detailed models of deuteron structure. 

Acknowledgements 

We are extremely grateful to J.L. Friar for insistently attracting our attention to the role of the 
3 Si positive-energy state in the inelastic contribution. The investigation was supported in part 
by the Russian Foundation for Basic Research through Grants Nos. 03-02-17612, 03-02-16510, 
and through Grant for Leading Scientific Schools. 



References 

[1] D.J. Wineland and N.F. Ramsey, Phys. Rev. A 5 (1972) 821. 
[2] G.T. Bodwin and D.R. Yennie, Phys. Rev. D 37 (1988) 498. 

[3] N.F. Ramsey, in Quantum Electrodynamics, ed. T. Kinoshita, World Scientific, 1990. 

[4] A. Bohr, Phys. Rev. 73 (1948) 1109. 

[5] F.E. Low, Phys. Rev. 77 (1950) 361. 

[6] F.E. Low and E.E. Salpeter, Phys. Rev. 83 (1951) 478. 

[7] D.A. Greenberg and H.M. Foley, Phys. Rev. 120 (1960) 1684. 

[8] A.I. Milstein, LB. Khriplovich, and S.S. Petrosyan, Zh. Eksp. Teor. Fiz. 109 (1996) 1146 
[Sov. Phys. JETP 82 (1996) 616]; Phys. Lett. B 366 (1966) 13. 

[9] F.E. Low, Phys. Rev. 96 (1954) 1428. 

[10] M. Gell-Mann and M.L. Goldberger, Phys. Rev. 96 (1954) 1433. 

[11] R.P. Feynman, Photon- Hadron Interactions (W.A.Benjamin, Inc., Reading, 
Massachusetts, 1972), p. 219. 

[12] A.I. Milstein, Siberian Journal of Physics (1995) No 1, 43. 

[13] A.C. Zemach, Phys. Rev. 104 (1956) 1771. 

[14] R.N. Faustov, A.P. Martynenko, |hep-ph/0211445 



